Structural arrest in dense star polymer solutions 
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The dynamics of star polymers has been investigated via extensive Molecular- and Brownian 
Dynamics simulations for a large range of functionality / and packing fraction r\. The calculated 
isodiffusivity curves display both minima and maxima as a function of f] and minima as a function of 
/. Simulation results are compared with theoretical predictions based on different approximations 
for the structure factor. In particular, the ideal glass transition line predicted by mode-coupling 
theory is shown to exactly track the isodiffusivity curves, offering a theoretical understanding for 
the observation of disordered arrested states in star polymer solutions. 

PACS numbers: 83.60.Hc, 64.70.Pf, 82.70.Dd, 83.80.Uv 



Star polymers play an important role in soft condensed 
matter systems, since they have been shown to interpo- 
late between hard colloids with a strong repulsive core 
on one side and the soft flexible polymeric systems on 
the other 0. Star polymers are constituted by a given 
number of polymeric arms, called functionality /, tied 
to a common center. As the functionality increases, the 
system becomes more similar to a hard sphere (HS) sys- 
tem, while lowering the functionality the inter-particle 
potential becomes increasingly soft. Recently, many in- 
teresting properties of star polymers have been clari- 
fied on the basis of an effective interaction potential be- 
tween star-polymer centers 0. In line with their hybrid 
polymer-hard colloid character, star polymers display no 
crystallization transition when the functionality / is low, 
/ < 34. At higher functionalities, a freezing transition 
takes place at about the overlap concentration of the sys- 
tem, into a bec-solid for lower functionalities and into an 
fec-solid for higher ones The freezing is succeeded 
by either a reentrant melting transition to the fluid for 
intermediate functionalities, 34 < / < 54, or by a cas- 
cade of structural phase transitions at higher values of /. 
The functionality-dependent bec- and fee- solids 0, 0|, as 
well as the reentrant melting transition have been 
experimentally observed in solutions of star-like block 
copolymer micelles. Though the crystalline solids are 
the phases of thermodynamic equilibrium at such high 
concentrations, the experimental situation is often some- 
what different. A variety of studies with star-polymers 
or star-like systems of various functionalities has shown 
that it is quite difficult to nucleate a crystal. Especially 
at high functionalities, the solutions display a gelation 
transition, i.e., a dynamical arrest into an amorphous 
crowded state in which the characteristic relaxation time 
of the system becomes extremely long [71 1^. Wi ITol ITU Xv\ . 
The onset of gelation, as opposed to crystallization, is 
further enhanced by the presence of some polydispersity 



in the samples, which gets indeed more pronounced as 
functionality increases. 

The purpose of this work is to analyze the dynamics of 
star polymers in athermal solvents theoretically, by em- 
ploying a combination of methods. Using the star-star 
effective interaction potential, in which all microscopi- 
cally fluctuating degrees of freedom are averaged out, we 
carry out detailed Molecular Dynamics (MD) and Brown- 
ian Dynamics (BD) simulations to measure the diffusivity 
of the star polymer fluids down to the homogeneous nu- 
cleation limit. Moreover . we carry out a Mode-Coupling 
Theory (MCT) analysis [Rj of the long-time limit of the 
correlation functions, which allows us to locate the non- 
ergodicity (ideal glass) transition line of the system. We 
find that, on increasing the number density of polymer, 
the characteristic time of the system goes through a se- 
quence of maxima and minima which we show to be re- 
lated to the oscillatory behavior of the effective HS diam- 
eter [l4|. We also find a strong correlation between the 
equilibrium phase diagrams of the system and the ideal 
glass line. The equilibrium reentrant melting transition 
is shown to have its counterpart in the reentrant melting 
of the disordered glass nested between two stable fluid 
phases. Finally we discover a striking similarity in the 
shape of the isodiffusivity and MCT ideal glass curves, 
regardless of the detailed approximation employed in the 
calculation of the structure factor S(q) and the type of 
dynamics (MD or BD). These findings strongly support 
the interpretation that the structural arrest of star poly- 
mers is a glass transition of 'effective hard spheres' char- 
acterized by a r\- and /-dependent HS diameter, despite 
the fact that the variety and nature of the equilibrium 
phases themselves is very different from the HS case. 

Standard MD and BD techniques have been employed 
in order to study the dynamic behavior of the star poly- 
mer system [l^. MD and BD lead to identical predic- 
tions for the long-time dynamics of the system within the 
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framework of the MCT lllUg. Thus > though BD offers 
a more realistic description of the short-time diffusion of 
the particles, the long-time behavior which is the relevant 
one for the glass transition is the same in both. Hydrody- 
namic interactions are ignored and the analysis is based 
on the effective center-center interaction potential 
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where (3 — 1/fcgT and r is the distance between the two 
centers. This potential is a combination of a logarithmic 
interaction at short distances, which gives the interac- 
tion its ultra soft character and stems from the scaling 
analysis of Witten and Pincus 01 1 an d a Yukawa form 
for the decay at long distances matched at a distance 
r = a. The quantity (3V(r) depends solely on a and 
/. The validity of this potential has been demonstrated 
via extensive comparisons both with small angle neutron 
scattering data [2J, [lOj and with monomer-resolved sim- 
ulations Jig. I n this entropic interaction, / _1 plays the 
role of T in normal fluids. 

In order to locate the line of structural arrest, the 
mean squared displacement averaged over all particles 
(|rj(f) — ri(0)| 2 ), where is the position of particle i, 
has been computed as a function of time t for several 
values of the two control parameters, the functionality / 
and the packing fraction 77 = irpa 3 /6, where p is the num- 
ber density. The values of the self-diffusion coefficient D 
are calculated from the relation 6Dt — (|rj(t) — r,(0)| 2 ), 
valid for large t. 

For the same system, the MCT equations 0] for mod- 
elling the structural arrest in supercooled liquid states 
have been solved; MCT furnishes a time evolution equa- 
tion for the normalized density fluctuation correlator 
Q(q,t) as a function of the momentum transfer q and 
t, which contains a term nonlinear in the correlator it- 
self. Knowing S(q) corresponding to the given interac- 
tion potential, the memory kernel entering in the non- 
linear term of the MCT equation can be evaluated and 
the equation solved for various values of q. In particular 
the non-ergodicity transition leading to structural arrest 
is obtained by performing the limiting value of the den- 
sity correlator lim^oo <fr(g, t) = /(g). The non-ergodicity 
factor f(q) is the solution of the equation 
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The vertex functions V, the coupling constants of the 




FIG. 1: The self diffusion coefficient, Z?md obtained by 
Molecular dynamics (upper panel, MD units E3)> £>bd by 
Brownian dynamics (middle panel, BD units) and Dus for 
an equivalent HS system (lower panel, arbitrary units), as 
a function of packing fraction for various values of the func- 
tionality [lyl . Missing points along the simulation curves were 
identified as crystalline states. 



theory, are 

V(q, k) = 4 [q • (q - k) c(|q - k|) + q • k c(k f x 
Q 

S(q)S(k)S(\q-k\) (4) 

and depend only on the Fourier transform of the direct 
correlation function c(q), or equivalently on S(q). In 
the A2 bifurcation scenario of MCT :ia the solutions 
of Eq. jump from zero to a hnite value at the ideal 
glass transition. The locus of the fluid-glass transition 
can be calculated varying the control parameters of the 
system, / and 77. 

S(q) has been calculated for the potential of Eq. us- 
ing two different approximate closures for the Ornstein- 
Zernike equation. The first one is the self-consistent 
Rogers- Young (RY) closure [2(| , which satisfies the ther- 
modynamic requirement that the isothermal compress- 
ibilities calculated either following the so-called fluctua- 
tion route or the virial route be identical to one another. 
The second approach is the modified hypernetted-chain 
approximation (MHNC) [2l|. In the MHNCjOne approx- 
imates the bridge function of the system [22j with that 
of an effective HS fluid. The optimum value of the di- 
ameter of this reference HS system is chosen in such a 
way to satisfy the Lado criterion |2^. Both the RY and 
the MHNC closures have been shown to yield results that 
compare extremely well with simulations [bH |24| . 
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FIG. 2: Experimental iso diffusivity curves for two values of 
the MD-diffusion coefficient (D = 0.067 and D = 0.017) and 
the MCT fluid-glass lines (computed with RY and MHNC). 
The equilibrium phase diagram (from Ref. |3j) is also reported 
for comparison. The inset compares the RY-ideal MCT glass 
line with the isodiffusivity data for D = 0.017, after a shift 
along the f~ l axis. 



The diffusion coefficient D for MD and BD are re- 
ported in Fig. ^ in the upper and middle panels respec- 
tively for a large set of / and 77 values In simple 
liquids, D decreases monotonically on increasing 77. In 
the present system, D has an highly non monotonic be- 
havior. At / = 32 a minimum around r\ = 0.5 is fol- 
lowed by a flat maximum up to 77 = 1.5. At larger / 
values (/ > 40) a clear sequence of minima and max- 
ima is observed. The MD and BD results show the same 
trend, and in supercooled states the D values become 
proportional. Such a proportionality supports the in- 
terpretation of the present results in term of MCT for 
supercooled liquids. Indeed, one of the basic prediction 
of the theory is the independence of the slow dynamics 
from the microscopic dynamics [T3I fljil |25| . 

The 77 and / dependence of D can be rationalized by 
calculating the isodiffusivity curves [2(| , i.e. the locus of 
points in the (77, /) plane where D has a constant value. 
Such curves, for two different values of D, are shown in 
Fig. [3 together with the equilibrium fluid-crystal coexis- 
tence lines, as calculated in Ref. The isodiffusivity 
curves run parallel to the coexistence lines, suggesting 
that crystallization in this system is achieved at the same 
D value, independently from the underlying crystalline 
phase, and offering thereby a strong confirmation of the 
dynamical freezing criterion of Lowen et al. j2?J • Fig. [21 
also shows the MCT ideal glass transition line (which 
can be considered as the isodiffusivity curve in the limit 
D = 0) using as input the RY and the MHNC S(q). Both 
the RY and the MHNC ideal glass lines track the molecu- 
lar dynamics isodiffusivity curves. The simulation curve 
can be perfectly superimposed to the theoretical curve 




FIG. 3: The MHNC structure factor S(q) (dashed line) and 
the non ergodicity factor f(q) (solid line) for three represen- 
tative points along the fluid-glass transition line. 



after a shift in as shown in the inset of Fig. [3 The 
agreement between theory and simulation below 77 = 1 is 
extremely good, independently of the / value. Even the 
reentrant shape of the isodiffusivity curve for / > 500 is 
captured by the theory. 

The shape of the ideal glass line in the (77, /) plane 
suggests interesting possibilities for modulating the dy- 
namics in star polymer systems. If crystallization can 
be avoided, a glass state can be generated by compres- 
sion which can then be melted with further compression 
(for example for / = 36) . Similarly, the dynamics of the 
polymer can be slowed down and increased again upon 
increase of its functionality (for example at 77 = 0.23). 

In order to better grasp the peculiar shape of the dy- 
namical arrest curve, we test the hypothesis that the 
slowing down of the dynamics is controlled by the 77 
and / dependence of an effective hard core. MHNC 
provides a well defined way for calculating the equiva- 
lent HS diameter cths and thus ?7hs- In the case of the 
star polymer potential, employing the Lado criterion |23| 
for the MHNC, the dependence of 7/hs on 77 has been 
already studied and it has been shown that it reflects 
the features of the interparticle interaction On the 
other hand, values of the diffusion coefficient Dhs for HS 
systems as a function of 77ns , the HS packing fraction, 
can be easily obtained directly by MD simulation. In 
a similar way -Dhs is known to follow for BD the sim- 
pie lawH MM Dm = A>(1 - 2t7 HS ) for 77ns < 0.4, 
with Do a constant. This offers the possibility of convert- 
ing the known 77-dependence of the effective HS packing 
fraction, 77hs(t7), into an effective -Dhs(?7) °f the corre- 
sponding HS system (both for BD and MD). This, in 
turn, can be compared with the star-polymer D values 
reported here. The comparison for MD is shown in the 
lower panel of Fig. ^ which illustrates that the simple 
mapping between the packing fractions of the star poly- 
mer 77 and the HS system 77ns captures the main features 
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of the diffusion coefficient, e.g., the location of the min- 
ima and maxima of the curves. The agreement between 
the two sets of data suggests that the slow dynamics in 
star polymer systems can be traced back - via a density- 
and functionality-dependent effective HS diameter - to 
the slow dynamics of the HS system (which is accurately 
described by MCTH3). 

Finally, Fig. [3| shows the non-ergodicity factor f(q) 
at three different representative points along the ideal 
glass line, calculated with the MHNC S{q). The f(q) 
shape changes continuously along the ideal-glass line, go- 
ing from the typical HS shape at small 77 (left and central 
panel) to a much more structured shape at large r\ (right 
panel). The /(g) width is a measure of the inverse of the 
cage localization length, which decreases on increasing 77. 
The q dependence of f{q) is always in phase with S{q), 
a feature common to all previously studied models. Two 
interesting features which develop at large 77 are the sig- 
nificant increase in f(q) in the g-rangc at qa = 7.44 and 
the small value for vanishing q. 

The MD and BD simulation data reported in this Let- 



ter, for a large range of / and 77, show that the dynamics 
of star polymer solutions is extremely rich. The isodif- 
fusivity curves have been shown to display both minima 
and maxima as a function of 77 and minima as a function 
of / which have been successfully connected to the be- 
havior to the 77 and /-dependence of the effective hard 
core diameter of an equivalent HS system. The detailed 
comparison between theoretical predictions and simula- 
tion confirms that MCT is a valid approach for guiding 
the interpretation of the disordered arrested states of soft 
matter materials j33, HH, offering a theoretical under- 
standing for the observation of disordered arrested states 
not only in colloidal systems characterized by a hard core 
[24l HE HE or charge-stabilized colloidal dispersions 
|38j but also in ultrasoft systems like star polymer solu- 
tions. 
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